Effective charge saturation in colloidal suspensions 
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Because micro-ions accumulate around highly charged colloidal particles in electrolyte solutions, 
the relevant parameter to compute their interactions is not the bare charge, but an effective (or 
renormalized) quantity, whose value is sensitive to the geometry of the colloid, the temperature or 
the presence of added-salt. This non-linear screening effect is a central feature in the field of colloidal 
suspensions or polyelectrolyte solutions. We propose a simple method to predict effective charges 
of highly charged macro-ions, that is reliable for mono-valent electrolytes (and counter-ions) in the 
colloidal limit (large size compared to both screening length and Bjerrum length). Taking reference 
to the non linear Poisson-Boltzmann theory, the method is successfully tested against the geometry 
of the macro- ions, the possible confinement in a Wigner-Seitz cell, and the presence of added salt. 
Moreover, our results are corroborated by various experimental measures reported in the literature. 
This approach provides a useful route to incorporate the non-linear effects of charge renormalization 
within a linear theory for systems where electrostatic interactions play an important role. 



I. INTRODUCTION 



When a solid-like object, say a colloidal particle (polyion), which carries a large number of ionisable groups at 
the surface is immersed in a polarizable medium (with a dielectric constant e, say water), the ionisable groups 
dissociate, leaving counter-ions in the solutions and opposite charges at the surface. The interactions between the 
charged colloids, which determine the phase and structural behaviour of the suspension, is mediated by the presence 
of micro-ions clouds. The complete description of the system is thus a formidable task in general. However in view of 
the large asymmetry of size and charge between macro- and micro- ions, one expects to be able to integrate out the 
micro-ions degrees of freedom, and obtain an effective description involving macro-ions only. In the pioneering work of 
Derjaguin, Landau, Verwey and Overbeek [1], micro-ions clouds are treated at the mean-field Poisson-Boltzmann (PB) 
level, yielding the foundations of the prominent DLVO theory for the stability of lyophobic colloids. An important 
prediction of the theory is the effective interaction pair potential between two spherical colloids of radii a in a solvent 
which, within a linearization approximation, takes the Yukawa or Debye-Hiickel (DH) form: 
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where Z is the charge of the object in units of the elementary charge e and n denotes the inverse Debye screening 
length. The latter is defined in terms of the micro-ions densities {p a } (with valences {z a }) as: k 2 — 4n£ B J2 a p a z 2 . 
The Bjerrum length £g is defined as £b — e 2 /(AirelcsT), where e is the permittivity of the solvent considered as a 
dielectric continuum: £b — 7 A for water at room temperature. 

However, this approach becomes inadequate to describe highly charged objects for which the electrostatic energy 
of a micro-ion near the colloid surface largely exceeds ksT, the thermal energy, because the linearization of the 
PB equations is a priori not justified. In this case however, the electrostatic potential in exact [2,3] or mean-field 
[4,5] theories still takes the Debye-Hiickel like form far from the charged bodies, provided that the bare charge Z 
is replaced by an effective or renormalized quantity Z c ^. The micro-ions which suffer a high electrostatic coupling 
with the colloid accordingly accumulate in its immediate vicinity so that the decorated object, colloid plus captive 
counter-ions, may be considered as a single entity which carries an effective charge Z c ff, much lower (in absolute 
value) than the structural one. Within the prominent mean-field PB theory [6] -often quite successful despite of its 
limitations-, Z and Z c s coincide for low values of the structural charge, but Z c g eventually reaches a saturation value 
Z^f independent of Z when the bare charge increases [5,7]. Arguably, the difference Z — Z c s is identified with the 
amount of counter-ions "captured" or "condensed" [8] onto the macro-ion. 

A reminiscent effect has been recognized in the physics of polyelectrolytes under the name of Manning-Oosawa 
condensation. Here, the object is an infinitely long and thin rod bearing A charges per unit length. At infinite 
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dilution and in the absence of salt, it can be shown at the PB level that the polyclectrolyte is electrostatically 
equivalent to a rod carrying A cqu ; v charge per unit length, where the equivalent charge density saturates to a critical 
value A cqu i v = 1/£b when A increases [9-11]. In general however, PB theory can be solved analytically in very few 
geometries only and the difficulty remains to predict Z e g for a given colloidal system [4,5,7,12,13]. In the absence 
of a general analytical framework for the computation of the effective charge, this quantity is often considered as an 
adjustable parameter to fit experimental data [14,15]. 

The aim of the present paper is to propose a method that allows to compute effective charges comparing favorably 
with PB in the saturation regime, provided the size a of the charged macro-ion is much larger than Bjerrum length 
£b and screening length k _1 . In the infinite dilution limit, we define the effective charges from the large distance 
behaviour of the electrostatic potential created by the (isolated) macro-ion [16]. While other definitions have been put 
forward [4,17-19] this choice is relevant in view of computing a macro-ion pair potential at large distances, in the spirit 
of the DLVO approach [20]. It moreover avoids the ambiguity of introducing a cutoff region in space which interior 
containing the colloid would exactly enclose a total charge equal to the effective one. At leading order in curvature 
(«a) _1 , our method easily provides effective charges at saturation close to their counterparts obtained in PB theory. 
In the situation of finite colloid concentration where it is no longer obvious to extract an effective charge from the large 
distance behaviour of the electrostatic potential computed within a non-linear theory, we follow the proposition put 
forward by Alexander et al. [5] introducing a Wigner-Seitz cell. In this situation, we generalize our original method 
into a prescription that we successfully test in various geometries, for different thermodynamic conditions (isolated 
systems or in contact with a salt reservoir). 

The paper is organized as follows. We first recall the basic framework of PB theory in section II. We then examine 
in some details the simple case of a spherical polyion in the infinite dilution limit (section III). This example allows 
us to devise a general method to compute the effective charge for arbitrary colloidal systems. The situation of finite 
density of colloids is then examined introducing Wigner-Seitz cells. The salt free case is developed in section IV, 
while the situation of finite ionic force is explicited in section V. We finally confront the results obtained within our 
prescription with experimental or simulation data in various geometries in section VI. We discuss the general validity 
of our mean-field treatment relying on PB approximation in section VII and conclusions are drawn in section VIII. 
The preliminary results of this study have been published elsewhere [21]. 



II. GENERAL FRAMEWORK: POISSON-BOLTZMANN THEORY 



Poisson-Boltzmann theory provides a mean field description of the micro- ions clouds in the presence of the polyions, 
acting as an external potential. The key approximation in the approach is the neglect of (micro-)ionic correlations. 
The size of the micro- ions with density p is neglected as well and the chemical potential reduces to its ideal contribution 
p = fcsTln(pA 3 ), where A is an irrelevant length scale. Without loss of generality, the macro-ions are supposed to be 
positively charged. 

At equilibrium the electro- chemical potential of the micro-ions is uniform over the system. Introducing the reduced 
electrostatic potential <f> = eV/ksT, the equilibrium condition for micro-ions thus reads at the mean field level 

ln(p ± A 3 ) ±<j> = ln(p A 3 ) (2) 

where are the density fields of the charged micro-species (counter-ions and co-ions), which we assume for 

simplicity mono-valent. The constant po will be specified hereafter. We restrict here to mono-valent micro-ions (both 
counter- ions and salt). For higher valences, the reliability of PB deteriorates (see section VII). The equilibrium 
condition, Eq. (2), is closed by Poisson's equation for the electrostatic potential 

V 2 = -Aw£ B (p+ - p~) . (3) 

The gradient of equation (2) expresses the condition of mechanical equilibrium for the fluid of micro- ions [22]. At this 
level, one has to separate between the no-salt and finite ionic strength cases. 

• No-Salt Case Only the released (here negative) counter-ions are present in the system. The PB equation for 
the reduced potential thus reads 

V 2 (f> = K 2 e* (4) 

where the screening constant k is defined as n 2 — Antspo, with p the constant introduced in Eq. (2). The 
latter is fixed by the electroneutrality condition, which imposes 
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with N c the number of (identical) macro- ions, of charge Ze, contained in the volume V. The quantity po is a 
Lagrange multiplier associated with the electroneutrality condition and has no specific physical meaning; it is 
modified by a shift of potential, which can be chosen at our convenience to fix cf) at a given point in the solution. 

• Finite Ionic strength situation In the finite ionic strength case, salt is added to the solution, so that both co- and 
counter- ions are present in the system. In the following we shall work in the semi-grand ensemble, where the 
{ colloids+micro-ions } system is put in contact with a reservoir fixing the chemical potential of the micro-ions 
po. In this case po in Eq. (2) is the concentration of salt in the reservoir (where </> is conveniently chosen to 
vanish), so that po = kBT\n(p A 3 ). Since we are considering mono-valent micro-ions po coincides with the ionic 
strength of the reservoir which is generally defined as Io = n~ 1 J2 a z aPa f° r a number n a of micro-ions species 
with valences z a and reservoir densities p° a . This results in the PB equation for the reduced potential <f> 

V 2 = K 2 sinh0 (6) 



where the screening factor k is now defined in terms of the micro-ion concentration in the reservoir k 2 = 
8n£ B Io = k 2 cs . 

In addition to these two situations, we shall also consider the case of infinite dilution where an isolated macro-ion 
is immersed in an electrolyte of given bulk salt concentration Io which thus plays the role of a reservoir. 

PB equations, (4) or (6), are supplemented by a set of boundary conditions on the colloids, expressing the relation- 
ship between the local electric field and the bare surface charges of the colloidal particles, ae. This gives the boundary 
condition for 4> at the surface of the colloid in the form 

(V0) • n = -4tt£ b o-, (7) 

where n denotes a unit vector normal to the colloid's surface. Except in simple isotropic geometries [23], the analytical 
solution of PB theory is not known. 



III. INFINITE DILUTION LIMIT: ASYMPTOTIC MATCHING FOR THE EFFECTIVE CHARGE 



In this section, after recalling a few results on the planar case, we explicit our method on the particular example 
of spheroids. We then generalize it to an arbitrary colloidal object and consider the case of charged rods as an 
application. We work in the infinite dilution limit, and therefore, we reject the external boundaries of the system at 
infinity. 



A. Planar Case 



In the case of the planar geometry, the non linear PB equation can be analytically solved. The detailed solution 
is given in appendix A. The important result however is that far from the charged plane, the solution of the PB 
equation reduces to that of the LPB equation : 



PB 



(z) ~ <p s e- KZ (8) 



The apparent potential (ps is equal to (ps — 4 in the limit of high bare charge of the plane. 

In this limit, the fixed charge boundary condition is therefore replaced on the plane by an effective fixed surface 
potential boundary condition 4> LPB (z = 0) = (f>s = 4. The effective charge density (in the saturation -high bare 
charge- limit) is then computed using Gauss theorem at the surface, yielding 

= (9) 
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B. Charged spheres 



Let us now consider a highly charged isolated sphere (bare charge Ze, radius a) immersed in a symmetric 1:1 
electrolyte of bulk ionic strength Iq. Within PB theory, the dimensionless electrostatic potential obeys equation (6). 
Suppose we know the exact solution (f> pB (r) (in spherical coordinates with the origin at the center of the sphere), 
and the bare charge Z is large enough so that the reduced electrostatic potential at contact, (f> pB (a), is (much) larger 
than 1. Then, we can divide the space surrounding the polyion into two sub- regions: a non- linear region (close to 
the particle's boundary) where </> pb (t") > 1, and a linear region where </> PB (r) < 1 (the potential vanishes at infinity). 
The surface delimiting these two regions is a sphere of radius r* such that PB (r*) — 1. 

Far from colloid, the complicated non-linear effects have died out to a substantial degree, and the solution also 
obeys the linearized Poisson-Boltzmann (LPB) equation V 2 </> = k 2 </>, and therefore takes the Yukawa form 

Z ff p-K(r-o) 

<W r ) = n • (10) 

LPB (1 + na) r 

The effective charge Z c g is defined here without ambiguity from the far field behaviour of <j> pB (r) : 

^^lpbW/^pbW = !• (11) 

In practice, LPB and 4> pB coincide in the linear region (r > r*), so that LPB ( r *) — 1 (i- e - is a quantity of order 
one). 

When a 3> k^ 1 , the non- linear effects are confined to the immediate vicinity of the macro- ion, with an extension 
kT 1 . We therefore have r*/a ~ 1 and as a consequence, LPB ( r = a) — ^lpb^ 7 "*) ~ 1" ^ c ^ us °btain the 
effective boundary condition that LPB is a quantity C of order one for r = a; from Eq. (10) this means that 
Z^g = Ca (1 + na) /£ B - This simple argument provides the non trivial dependence of the effective charge at saturation 
upon physico-chemical parameters; it applies in the saturation regime of PB theory where Z c g — Z^g and assumes 
that the bare charge Z is high enough so that the non- linear region exists. In order to determine the constant C, wc 
may consider the planar limit a — > oo where the analytical solution of PB theory is known (see above and appendix 
A): the surface charge density Z^ / (Ana 2 ) should coincide with that of a charged plane k/(tt£ b ), Eq. (9). This 
imposes that C = 4 and going back to the charge: 

An 

Z%=T (l + «a). (12) 

In de-ionized solutions, this argument leads to the scaling Z^ cx a/is, which has been recently tested for various 
latex colloids [24]. 

The physical argument leading to (12) may be rationalized as follows. The situation of large na corresponds to a 
low curvature limit where the solution of Eq. (6) may be approximated by the solution of the planar problem in the 
region where curvature effects may be neglected : the latter corresponds to a region a < r < a + 8a, with 8a ~ a. 
It is crucial to note that r* < a + 8a since, as mentioned above, the extension of the region where the non linear 
effects are important (defining r*) has an extension of order k" 1 , smaller than 8a ~ a in the limit of large na. As 
a consequence, in the region r* < r < a + Sa, the solution of the LPB equation, Eq. (10), may be matched to the 
asymptotic expression of the planar solution, given in Eq. (8) (using r ~ a and z ~ r — a). The expression (12) is 
therefore recovered, showing again that at the linearized level, the apparent potential is LPB (a) = 4 in the saturation 
limit. 

Eq. (12) provides by construction the correct large na behaviour of Z^ , and becomes exact (compared to PB) 
in the planar limit. We will show below that it remains fairly accurate down to na of order 1. A similar expression 
may be found in [25,26], but the generality of the underlying method does not seem to have been recognized. This 
result is supported by the work of Oshima et al. [27] which proposes an approximation scheme of the non-linear PB 
equations for spheres in infinite dilution, for large na. In particular these authors obtain an analytical approximation 
for the apparent potential at the colloid surface, which reads in the saturation regime: 

^Osh =8 l + ^ (13) 

Supplemented with expression (10), this leads to the improved effective charge 

t _ 8a (1 + na) 2 

ZcS " T B \ + 2na ' (U) 
In the limit of large na where (j>s — > 4, both Eqs. (12) and (14) have the same behaviour. 
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FIG. 1. Effective charge in the saturation regime Z^Ib/cl as a function of Ka for spheres in the infinite dilution limit with 
added salt. The symbols (open circle) are the "exact" solution estimated from the large distance behaviour of the electrostatic 
potential solution of the full non-linear PB equation. The continuous (resp. dashed) line is Z sa t found with Eq. (12) [resp. Eq. 
(14)]- 



In order to test the validity of these results, we have numerically solved the full non- linear PB equation, Eq. (6) 
and computed the effective charge from the electrostatic potential at large distances, i.e. the value required to match 
^lpb to tne ^ ar ^pb b tame( i numerically. For each value of Ka, we make sure to consider large enough bare 
charges in order to probe the saturation regime of Z e g. Figure 1 compares the numerical PB saturation value of the 
effective charge to the prediction of our approach, Eq. (12), and to that obtained using the results of Oshima et al, 
Eq. (14). We see that Eq. (14) provides an accurate estimate for as a function of Ka, for na > 1. Working at 
the level of our approach only, Eq. (12) still yields a reasonable estimate for Z*g(na), specially for high values of the 
parameter na. In the limit of small Ka, both expressions (12) and (14) differ notably from the PB saturation charge 
which diverges, as shown by Ramanathan [28], as: 

Z eS ~ IT I -2 H Ka ) + 2 M- M«°)] + 4 In 2} for Ka -» 0. (15) 

At this point, it is instructive to briefly reconsider the work of Squires and Brenner [29] who demonstrated that the 
attractive interactions between like-charged colloidal spheres near a wall could be accounted for by a non-equilibrium 
hydrodynamic effect (see also [22]). In their analysis, they used an ad hoc value of 0.4 for the ratio Cgiass/csphcrc °f 
surface charge densities of planar and spherical polyions, in order to capture the one-wall experiment of Larsen and 
Grier [30]. This was the only free parameter in their approach. From Eqs. (A6) and (14) for the saturation values, 
we easily obtain 

ggiass _ Ka(l + 2ko) 

Csphere 2(1 + Ka) 2 

In the experiment of [30], we have Ka ~ 1.2 [larger than 1, so that (14) is reasonably accurate], and we obtain 
Ogiass/osphcrc — 0.42: it is thus possible to justify the choice made in [29] assuming that both the confining wall and 
the pair of colloids are charged enough to sit in the saturation regime. In this respect, knowledge of their bare charges 
is unnecessary. 



C. Arbitrary colloidal object 

Generalizing this analysis for an arbitrary colloidal object (of typical size a), we propose the following method to 
estimate the effective charge in the limit of large values of Ka: 

1. Solve the LPB equation for the geometry under consideration 

2. Define the saturation value, Z^g, such that the linear reduced potential at contact is a constant, C, of order 
unity 

|0s - 0buik| = C, (17) 
where the asymptotic matching with the planar case yields C = 4. 



5 



3. If one is interested in the effective charge for arbitrary and possibly small bare charges, a crude approximation 
follows from 



z eS = z z < 

Z eS =Z s e $ Z^Zlf. (18) 

Our approach has several advantages. First, we do not need to solve the full non-linear PB equations to obtain 
the effective charge. Second, the proposed method provides an analytical prediction for Z^g. Third, our approach is 
easily adapted to other macro- ion geometries or finite dilutions, unlike that of Ref. [27] (even if these authors could 
find an equivalent of expression (13) for cylinders, see below). 

In the following, we will mainly focus on the high bare charge limit of the colloids where the effective charge reaches 
a saturation plateau, Z^. In order to simplify notations, we will denote this saturation value Z sa t- 



D. Rod-like macro-ions 




FIG. 2. Effective line charge density, £sA sa t, versus na (the reduced Debye-Hiickel constant) for cylinders in the infinite 
dilution limit with added salt. The symbols (open diamonds) are computed from the large distance behaviour of the electrostatic 
potential solution of the full non-linear PB equation, solved numerically. The continuous (resp. dashed) line is our estimate 
for ^BAsat, Eq. (21) [resp. the improved estimate, Eq. (23)]. 



Now the object is an infinitely long cylinder (radius a, bare line charge density Ae). The solution of linear PB 
equation is (in cylindrical coordinates where r is the distance to the axis) 

^)=2A^^§M, (19) 
na Ki(Ka) 

where K and K x are (respectively) the zero and first order modified Bessel functions of the second kind. Hence the 
apparent potential is 

4>s = 2M B ±^ (20) 

Setting (frs = C = 4 yields our estimate for the effective line charge density at saturation 

_ 2na KjjKa) 

Asat — —„ T7-7 r • (/J- J 

In the limit of large values of the bare line charge density A, Oshima et al. obtained an approximate expression 
for the apparent potential in the saturation regime (based on an approximation scheme for the PB equation, see 
Appendix A of Ref. [27]): 
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6 Osh = g Kjjna) 
S [K ( K a) + K^na)] 



(22) 



As expected, we note that (j)g sh — > 4 in the limit ko> 1. However, from Eq. (22), we deduce an improved estimate 
of A sat 



A, 



l sat 



i B Ko(ko) [K (ko) + K^na)]' 1 ' 

In Fig. 2, we display ^sA sa t [estimated either with Eq. (21) or Eq. (23)] as a function of na, together with the 
"exact" value of fsAsat found by solving the full non-linear PB equation for high bare charges in the saturation 
regime. Note that the plot is in log-linear scale, in order to emphasize the small na region where our method is not a 
priori expected to work. Surprisingly, the agreement between the numerical result of the full PB equation and (23) 
is satisfactory down to very low values of na, na ~ 10~ 2 , although the two quantities have a different asymptotic 
behaviour: the exact ^sA sat is finite when na — > (^sAsat = 2/7r from rcf. [31], see next paragraph), whereas both 
estimate Eq. (21) and (23) vanish, although extremely slowly (as — 1/ log(«;a)). 

Importantly, na — > is the asymptotic regime where the celebrated Manning limiting law [11,32] happens to be 
exact, and the condensation criterion holds. In this limit, above the condensation threshold, the electrostatic potential 
solution of the full non-linear PB equation is indistinguishable from that of a cylinder carrying a line charge density 
Acquiv = [11]- The two quantities A oqu j v and A sat may be coined as "effective charges", but we maintain our 

initial definition of the effective charge from the far field potential solution of the non-linear PB equation. In this 
respect, A cq uiv 7^ A sa t, (as already discussed in the Appendix of reference [11]). This is because one expects a remnant 
non-linear screening of A oqu i v , so that A sat < A oqu i v = 1/Ib- The limiting situation na — > has been solved recently 
within Poisson-Boltzmann theory, using exact results from the theory of integrable systems [31]. The corresponding 
solution illustrates our point. This seminal work allows to compute explicitly the effective charge, which reads 



2 . ( 7T 



lim A c ff = —r- sin 



in (^\l B ) . (24) 



Accordingly, when A exceeds the Manning threshold 1/Ib, the effective charge saturates to a value 

2 0.6366 1 

Asat — 7 — 7 < Acquiv — ~„ • (25) 

-KlB ZB t-B 

It is noteworthy that the limit 2/(tt£b) (compatible with the numerical results reported by Fixman, see for example 
Fig. 1 of rcf. [32]) is reached extremely slowly as na is decreased, in practice for na < 10~ 6 . For example, for 
na = 10~ 2 , the numerical solution of the PB equation yields A sa t = 0.81/^s, hence a value 30% larger than the 
asymptotic limit. 

IV. EFFECTIVE CHARGE AT FINITE CONCENTRATION. THE NO-SALT CASE 

The situation of finite density of colloids does not allow to define an effective charge from the far field of a single 
body potential, as done in section III. Here, we rely on the proposition put forward by Alexander et al. to define an 
effective charge [5]. We recall here the main points of this PB cell approach. First, the procedure makes use of the 
concept of Wigner-Seitz (WS) cells: the influence of the other colloids is accounted for by confining the macro-ion 
into a cell, with global electroneutrality [33-36]. The size of the cell, Rws is computed from the density of colloids, 
while its geometry is chosen as to mimic the spatial structure of the colloids in the solution. Second, the "effective" 
potential solution of the linearized PB equation is such that the linear and non-linear solutions match up to the 
second derivative at the boundaries of the WS cell (hence they match up to at least the third derivative because 
of electroneutrality in "isotropic" -spherical or cylindrical- cells). Note that in the original paper of Alexander et 
al. the procedure was introduced to obtain the effective charge from the numerical solution of the non linear PB 
equation. But in the present work we shall use the approach to get effective charges at the LPB level together with our 
prescription. Such a route has proven successful for mono-valent micro- ions, see e.g. [37-39], and it has been shown 
recently that similar ideas could be employed to describe discrete solvent effects (again for mono-valent micro-ions 
[40]). 

In this section, we generalize the analysis proposed in section III C to find a prescription suitable to treat the case 
of finite concentration of colloids. We eventually compare our results to those obtained following ref. [5], for planar, 
cylindrical and spherical geometries. 
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A. Generalized prescription and planar test case 

In the infinite dilution case, the reference potential is the bulk one ^buik- The natural generalization of this choice 
for the finite concentration case consists in replacing in Eq. (17) </>t>uik by <p s the reduced electrostatic potential at 
the boundary of the WS cell. Hence, we propose 

\<j>s-^\=C. (26) 

If added salt was present in the suspension (see section V), we should recover Eq. (17) from Eq. (26) in the infinite 
dilution limit where Rws goes to infinity. We consequently expect the value C = 4 to be relevant for the situation of 
finite density of colloids with added salt. Searching for a unified description, we also test the possible validity of the 
choice C = 4 in the no salt situation. It is therefore instructive, as an illustration of the method and benchmark, to 
analyze the simple case of a charged plane confined in a WS cell, without added electrolyte. As recalled in Appendix 
B, the analytical solution of the non linear PB equation is known in such a geometry when counter- ions are the only 
micro- ions present, which allows to check the validity of our assumptions in the limiting case of finite concentration. 
Below, we compare these "exact" results to the predictions of our prescription. 

The exact apparent potential, <fis is obtained using Eq. (B5) at x = for the plane: (f>s = cosh( -K" LPB h) — 1. 
Our prescription imposes (f>s = 4, yielding K hpB h — ArcCosh(5). The effective charge is obtained from Gauss 
theorem at the surface, i.e. Eq. (B6) with a replaced by cr e ff. This leads to the final result of our prescription 
<7 sat = \/6ArcCosli(5)<7c = 5.6cr c (where a c = l/wlsh), which should be compared to the exact result a sat — 5.06cr c [see 
Eq. (B12) in appendix B]. First it is striking to note that our prescription predicts the correct functional dependence 
of the effective charge in terms of the parameters of the system. Moreover the numerical prefactor in front of a c is 
only within 10% of the "exact" value obtained in Eq. (B12), which is quite a satisfactory agreement. 

However certainly the most interesting feature which comes out from the previous results is the fact that the 
apparent potential at contact, 4>s, obtained within the analytical resolution of the PB equation, does saturate to a 
constant value 4>s — 3.66 in the limit of very large bare charges: this value is very close to the value we prescribe, 
4>s = 4! This is a non trivial point, since the physical conditions in the present case are very different from the isolated 
plane case (previous section). We conclude that the analytic results available for a confined one-dimensional electric 
double-layer support our prescription. For a more refined analysis of the electrostatics of counter-ions between planar 
charged walls, going beyond PB, we refer to the work of Netz et al. [41]. 

In the remaining of this section, we further test our prescription against results for spherical and cylindrical macro- 
ions. 



B. Spheroids 

Here, the object is a charged spherical colloid (bare charge Ze, radius a) confined with its counter-ions in a concentric 
WS sphere (radius Rws)- The packing fraction is defined as rj = (a /Rws) 3 - PB equation is again linearized around 
the boundary of the WS cell, yielding Eq. (B4) which we recall here: 

V 2 = ^ pB (0+i). (27) 

As for the planar case, the boundary conditions are V(j)(Rws) = (clcctroneutrality) , LpB (i?) = (because we 
impose by commodity the potential to vanish at the WS cell, see Appendix B). The solution <fi LPB thus reads: 

e lpb e lpb 

W»0 = -1 + /+^— + f-—T- (28) 

with 



K, Rw s i 1 

f± = LP 9 \ eMTK hpB Rws)- (29) 

ZiV LPB 



The charge Z e g of the colloid is obtained from the spatial derivative of at the colloid surface: 

Z eS = ^--n- 1(1 - K* pB aRws)smh[K LpB (R W s - a)} - K hpB (R WS - a) cosh[K hpB (Ri W s - a)}\. (30) 
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At this level the screening constant K hpB is still unknown: it is fixed by our prescription which imposes the apparent 
potential of the colloid, such that 4 l L p B ( r = a) = C = A with 4> LPB ( r ) given in Eq. (28). The effective charge, Z sat , 
is eventually computed using Eq. (30). On the other hand K pQ is defined from the PB counter-ion density at WS 
boundary p~(WS): 

K 2 pB =^e B p-(WS). (31) 

The equation for K hpB , (f> LPB (r = a) = 4, is solved numerically using a simple Newton procedure. Fig. 3 displays 
the corresponding Z sat as a function of r\ together with the effective charge (again at saturation) found by solving 
numerically the full non linear PB equation, together with Alexander's procedure. We recall that this procedure 
defines the effective charge entering LPB equation such that the solution of PB and LPB equations match up to the 
second derivative at the WS boundary. We emphasize that once PB equation has been solved numerically, no further 
numerical fitting procedure is required to match LPB to (f> pB and compute the effective charge: the counter-ion 
density at WS boundary p~(WS) is known and K pB follows from Eq. (31). Replacing K LpB with this value in Eq. 
(30) then gives the "Alexander" Z c g (a similar remark holds with added salt, see below). We see again that our 
prescription works reasonably well. 
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FIG. 3. Effective charge at saturation versus packing fraction for spherical polyions without added salt. The symbols (open 
diamonds) represent the effective charge found by solving numerically the non linear PB theory supplemented with Alexander's 
procedure [5] . The continuous line is Z sa t within our prescription. 



C. Cylinders 

Here we apply the previous procedure to an infinite cylinder (radius a, bare charge per unit length Ae) enclosed in 
a WS cylinder (same axis, radius Rws)- We define the packing fraction as r\ = (a/Rws) 2 - This case is particularly 
interesting since the solution of the PB equation in the no-salt case is known from the work of Fuoss et al. and Alfrcy 
et al. [10]. This therefore provides another critical test of our prescription. We note that a similar approach has been 
followed by H. Lowen [38]. 

The calculation follows the same lines as for the previous spherical case. The solution of LPB equation (27) in the 
two dimensional case, with the usual boundary conditions and the choice </> LPB (-Rivs) = reads: 

^lpf» = -^ + K LPB R WS {h(K hpB R ws )K (K LpBP )+K 1 (K LpB R ws )I (K LpB p)}, (32) 

where use was made of the identity x[Iq{x)K 1 (x) + Ii(x)K (x)} = 1. From the spatial derivative of 4> LPB at r — a we 
deduce the effective line charge density A e fl=: 

Wb = ^Kl pB aR ws ^l 1 (K LpB R ws )K 1 (K LpB a) - I^K^a) Kr^^s)}- (33) 

If the above expression is evaluated replacing K LpB by the exact K pB following from (31) once the non- linear 
problem has been numerically solved, we obtain the original Alexander value (not necessarily at saturation, and 
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without having to implement in practice a numerical fitting procedure). On the other hand, as in the spherical case, 
the screening constant K LPB at saturation is obtained (approximately) by imposing the potential at the polyion's 
surface: </> LPB (a) = C = 4 in the previous equation for 4> LPB ■ Evaluation of Eq. (33) then gives the saturation value 
A sa t- 




0.1 0.2 0.3 0.4 



FIG. 4. Same as Fig 3 for charged rods, except that the non linear PB results are analytical here [10]. 

This result is compared with the effective charge deduced by applying Alexander's procedure to the analytical 
results of Fuoss et al. and Alfrey et al. in the large bare charge (saturation) limit [10]: the corresponding "exact" 
value for the effective charge is chosen such that the solution of the linearized PB equation, Eq. (32), matches the 
solution of the non linear PB equation (Fuoss/Alfrey et al. solution) at the WS boundary up to the second derivative. 
The resulting effective charge is plotted on Fig. 4 together with the value of the effective charge obtained within our 
prescription. In Fig. 5, we compare the screening factor K hpB obtained within our prescription to the "exact" value 
K pB derived from the analytical solution of the PB equation and Eq. (31) (again in the limit of a large bare charge 
of the cylinder where the effective charge saturates). The agreement between both quantities is remarkable, even up 
to extremely high packing fractions (80% on Fig. 5). 
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FIG. 5. Comparison of the exact inverse screening length K pB obtained from the solution of PB equation (diamonds) with 
its counterpart K^ pB obtained within our prescription (continuous line), for highly charged rod- like polyions without added 
salt. The inset shows the ratio K LpB /K pB as a function of packing fraction. 

Another interesting check concerns the apparent potential at the surface of the cylinder. Applying again Alexander's 
procedure to the exact solution of Fuoss/Alfrey et al. , one obtains the LPB potential which matches the exact solution 
up to its third derivative at the WS cell boundary. The value of this potential at the surface should be compared 
with the value we prescribe, i.e. (f>s = 4. The result is plotted on Fig. 6, showing again a good agreement except 
at low volume fraction, as expected (since as discussed in section III, our prescription is not expected to work in the 
very small na limit). 
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FIG. 6. Dependence on volume fraction of the reduced linearized contact potential, <f>g — 4>(a), with <f> the LPB potential 
matching the analytical solution of the PB equation, following Alexander's procedure. 



Finally, we report an intriguing result: in the limit of vanishing density it can be shown analytically that K m Rws 



\/2 [10]. Using this result together with K 



LPB 



0, we obtain from Eq. (33) 



lim A sat 



0.6358—. 



(34) 



This asymptotic value is displayed in Fig. 4 with an arrow. We observe that this limit is approached (although very 
slowly) as r\ decreases. Surprisingly, the result of Eq. (34) is very close to the exact expression (25) of Tracy and 
Widom [31] where the limit na —* is taken after that of infinite dilution. In principle, the limits of infinite dilution 
and of vanishing added salt have no reason to commute. The difference between the two A sat quantities illustrates 
this point, with the surprise that the results are nevertheless very close numerically: 



lim 

no salt 



lim A sat 

00 dilution 



lim lim A sat 

00 dilution no salt 



0.6358 — 
0.6366 -5-. 

*-B 



(35) 
(36) 



V. EFFECTIVE CHARGE AT FINITE CONCENTRATION. THE FINITE IONIC STRENGTH CASE 

We now turn to the case where salt is added to the colloidal suspension. More precisely, as already discussed in 
section II, we consider the semi grand canonical situation where the colloidal suspension is put in contact with a 
reservoir of salt, through a semi-permeable membrane (dialysis experiment). The concentration of mono-valent salt 
micro-ions in the reservoir po fixes the chemical potential of the micro-ions in the suspension. However, due to the 
presence of the charged macro- ions, the salt concentration in the solution, p Sl differs from that in the reservoir p : 
this is the so-called "salt exclusion" or "Donnan effect" [36,42,43]. 

As in the previous section the effect of finite concentration is accounted for within PB cell theory (using a WS 
sphere of radius Rws)- Here again, we use the prescription Eq. (26) to predict the effective charge of the macro-ions. 
For this purpose, it is convenient to choose that the electrostatic potential cf) vanishes in the reservoir so that PB 
equation reads 

V 2 = K 2 cs sinh</> (37) 

where the screening factor n rcs is defined in terms of the ionic strength of the reservoir: k 2 0S = 8itIbIo- 

Let us now consider the linearized ( "LPB" ) version of this equation. We again linearize around the value of the 
potential at the boundary of the WS cell, <f> s = (p(Rws), often referred to as the "Donnan potential" 

V 2 6<j> = K 2 LpB (S^ + l0 ) (38) 

where we introduced S(f> — (f> — <p s and: 
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70 = tanh[</> s ] = Wl 



LPB 



(39) 
(40) 



The second order differential equation (38) is solved invoking the two self consistent boundary conditions 

5<t> = and ^ = for r = R ws , (41) 
or 

so that 8(f) is known as a function of distance and depends parametrically on K hpB . We emphasize that K LPB is still 
unknown at this point. It is computed as in section IV from our prescription on the reduced potential 



4. 



(42) 



A. Spheroids 



N 




« 20 



FIG. 7. Effective charge (at saturation) of spherical colloids (radius a) as a function of volume fraction r\ for K ros a = 2.6. 
The continuous line is the effective charge (at saturation) computed using the prescription, while the symbols are the results 
of the non-linear PB cell theory, following Ref. [5]. 



With the same notations as above, the appropriate solution of LPB equation (38) is 

'LPB (0=70 -1 + /+— + /- 



LPB 



(43) 



where the functions f± are defined in Eq. (29). Note that expression (28) is recovered by taking the formal limit 
Kros = in the previous equation. 

Our prescription allows to compute K LpB at saturation, such that <54>(a) = 4, without any reference to the solution 
of the non linear PB problem. This equation is solved numerically using a Newton procedure. Once K BpB is known, 
the effective charge follows from the gradient of 8<j)(r) in Eq. (43) taken at r — a (it may also be computed by 
integrating the corresponding LPB charge density over the volume accessible to the micro-ions, i.e. a < r < Rws) 



Z sa .t — 7o 7— T7 
Ib A. 



lpb 



(K^ pB aR W S - 1) sinh[if LpB {Rws ~ a)} + K LpB (Riws - a) cosh[if LpB (Riws - a)} \ (44) 



1 — {K-rcs/K, PR ) 4 . Again, our prescription 84>s 



4 provides a value for K LpB which is 



with 70 = 

an approximation for the exact K pB at saturation, related to microions densities at the WS boundary through the 
expected Debye-like form 



Kb = 47:1 B [P + (WS)+ P -(WS)]. 



(45) 
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If Eq. (44) is evaluated with the exact K pB , Alexander's original effective charge follows (hence without having to 
implement any numerical fitting procedure). We also emphasize that as in the previous sections, the right hand side of 
Eq. (44) provides the effective Z e g a la Alexander (i.e not necessarily at saturation), once evaluated with the correct 
K pQ (deduced from the numerical solution of the non- linear problem). 

e-o r- nq . . =r- 

5.5 - 

5.0 - 

-©-4 5- D 
to ^ c 

4.0 



1.0 O o O 



0.9 1 1 1 

0.1 0.2 0.3 



3.5 - 



0.1 0.2 0.3 

FIG. 8. Dependence on volume fraction of the "exact" reduced linearized contact potential, 5<j>s = 4>{a) — (j>(Rws), with 
<j> the LPB potential matched to the numerical PB solution according to Alexander's procedure. Our prescription assumes a 
constant value, 5<f>s = 4. Inset: ratio K LpB /K pB versus packing fraction. The ratio is seen to be very close to unity over the 
explored packing fraction window. 



The results for the effective charge at saturation Z sat as a function of volume fraction are displayed on Fig. 7 for 
K Ics a — 2.6. As in the previous sections, we compare this result with its Alexander's counterpart at saturation. Our 
predictions are seen to be compatible with those obtained in the PB cell model. 
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FIG. 9. Influence of the point r match chosen to match the analytical LPB and numerical PB solutions on the effective charge 
in the saturation regime. The situation is that of a spherical polyion in a spherical WS cell, at packing fraction r\ — 0.05 and 
Kres a, = 2.6. 



From the numerical solution of the PB equation, it is possible to extract the apparent surface potential, 5<fts = <M r = 
a) — (fi(Rws) (in the latter expression (f> is defined as the solution of the LPB equation matching the full -numerical- 
PB equation up to second derivative at the WS cell boundary). By construction, this potential may be obtained 
inserting the numerically obtained K pB = k ics cosh 1 ^ 2 into (43). This apparent potential should be compared 
against our prescription 8<j) = 4. The corresponding result is shown on Fig. 8. We observe that 5(j)s indeed saturates 
to a value close to 4. The inset shows K hpB / K pB versus i], where K pB is the "exact" screening length for the LPB 
equation at saturation, obtained numerically; K is the same quantity estimated from our prescription. We observe 
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that although for small packing fractions 5<fis slightly departs from our approximation C = 4, the estimated K LpB is 
still remarkably close to the exact one. 

Independently of our prescription, we finally test the relevance of Alexander's procedure [5] in the following way. Z sat 
has been obtained above from the matching of a generic LPB potential to the numerical PB one at r matc h = Rws ■ It is 
also possible to implement the matching at a different location inside the cell, and we denote Z S at(f match) the associated 
effective charge, at saturation. This quantity, normalized by the "usual" one Z sat (Rws) is displayed in Fig. 9. For 
the packing fraction of 5% considered, r ma ,tch/Rws is bounded below by a/Rws — 0.37, and Z sa , t {Rws) — 16.7cl/£b, 
see Fig. 7. We observe that Z sat is relatively insensitive to r matc h for 0.6Rws < f match < Rws- 



B. Rod-like polyions 

Using the same notations as in section IV, the appropriate solution of the LPB equation in cylindrical geometry 
reads 



u>bW =lo<K hpB R ws 



h(K hpB R ws ) K„(K hpB p) + Ki (if LpB iVs) I (^ LPB p) 



(46) 



with 70 = ^Jl — (K res /if L p B ) 4 . Again, K LpB is obtained as the solution of the equation <5</> LPB (p = a) = 4. Once this 

equation is solved, the saturation value of the effective charge, A c ff , follows from the spatial derivative of the potential 
<5</> LPB at the rod surface: 

Acff - ^^LPB aiV ^° j I i(^LPB ffi ^s)K 1 (if LpB a) - I 1 (if LpB a)K 1 (^ LpB i?z ws )|. (47) 

As in the previous sections, Eq. (47) gives analytically the relation between the effective charge a la Alexander et al. 
and micro-ions densities at the WS boundary. As such, it applies for any value of the bare charge A, and in particular, 
for A — > 0, K^ pB is such that A e ff/A — > 1. A similar remark applies for Eqs. (30), (33) and (44). If we choose for 
K LpB the "exact" K pB value, we recover the "exact" cell model (Alexander) effective charges. However, the quantity 
K hpB solution of <50 LPB (p = a) = 4 is supposedly the inverse screening length at saturation and therefore provides 
an approximation of A sa t once inserted into (47). 

The corresponding results for A sa t as a function of volume fraction are displayed in Fig. 10 for n rcs a = 3. As in 
the previous sections, we compare this result with its counterpart obtained from the numerical solution of PB theory 
together with Alexander's procedure for the effective charge in the saturation limit. The agreement with the numerical 
results of the full non linear PB equation is seen to be satisfactory. 
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FIG. 10. Effective charge lsA sa t as a function of packing fraction for cylinders with added salt (K les a = 3). The symbols 
represent the effective charge at saturation within the usual PB cell approach, while the continuous line follows from our 
prescription. 
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VI. CONFRONTATION TO EXPERIMENTAL AND NUMERICAL RESULTS 



In the previous sections, we have tested our results for the effective charges against the numerical solutions of PB 
theory. However, the effective charge is a difficult quantity to measure directly in an experiment (see however the work 
reported in [24] confirming the scaling Z sat oc o/fe for low ionic strength suspensions of spherical latex colloids). In 
order to assess the experimental relevance of the above ideas, we now turn to the computation of osmotic properties 
for spherical and rod-like macro-ions, easily accessible both experimentally and within our approach. In the case of 
spherical colloids, we start by considering the phase behaviour of the suspension as a function of added salt. 

A. Crystallization of charged spheres 

The phase diagram of charged spherical colloids has been widely explored experimentally, in particular by 
Monovoukas and Gast [14]. In this work, the macro- ions were charged polystyrene spheres, with radius a ~ 660 A. 
The authors moreover compared their experimental phase diagram to that computed for particles interacting through 
a Yukawa potential (1) (the Yukawa phase diagram has indeed been investigated extensively by numerical simula- 
tions [44-46]). However such a comparison experiment/theory requires an ad-hoc choice for the effective charge Z e g 
[prefactor of Eq. (1)]. The authors found that a reasonable agreement with the numerical results was obtained for 
a specific choice of the effective charge, Z SSLt = 880 (although they reported conductimetry experiments indicating a 
macro-ion charge around 1200). 

We focus in the following on the melting line of the phase diagram obtained in [14] . We use here our predictions for 
the effective charges at saturation: we do not need to know the bare charge of the polystyrene spheres, as this quantity 
is presumably much larger than the corresponding saturation plateau of Z e g, which means, within the PB picture, 
that Z e fi ~ Z sat . Once Z sat (and the corresponding screening constant K LpB , see previous section) is known for a 
given density and ionic strength, it is possible to insert it into the computed generic phase diagram of Yukawa systems 
[44] to obtain the corresponding stable phase. We extract in particular the melting line from these numerical results: 
we prefer to use these numerical results for the phase diagram (instead of performing a full theoretical analysis) since 
our main focus remains to check the relevance of our predictions for the effective charges. This requires a "reliable" 
description for the melting line, which numerical simulations provide once the potential is given. 
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FIG. 11. Liquid-solid transition of charged polystyrene colloids: volume fraction for melting rj m as a function of salt ionic 
strength Jo- Dots are experimental points for the melting line extracted from Ref. [14]. The solid line is the theoretical 
prediction for the melting transition using our prescription for effective charges (see text) while the dashed line corresponds to 
an ad hoc fixed effective charge Z e g = 880, as proposed in [14]. 

We emphasize that at this level, the only parameter entering our description is the diameter of the beads, which 
is measured independently. Accordingly, there is no adjustable parameter in our equations and the resulting phase 
diagram is strongly constrained. The results for the melting line using our prescription for the effective charge are 
confronted to the experimental data in Fig. 11. We also plot the result for the melting line for an ad hoc constant 
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effective charge, Z sat — 880, as was proposed in Ref. [14]. The observed agreement supports the pertinence of our 
prescription for Z sat which reproduces the experimental phase diagram. In our case, the effective charge does vary 
between 500 and 2000 along the melting line, depending on ionic strength and density. This could explain that the 
conductimetry measurements performed independently by Monovoukas and Gast (although at an unspecified ionic 
strength) yield another value for the effective charge of the spheres (Z <~ 1200 as quoted above). 

B. Osmotic pressure of a suspension of spherical colloids 
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FIG. 12. Reduced osmotic pressure versus volume fraction for spherical polyions in the salt free case 

where lo (and thus K ros ) vanishes. The symbols are the PB values obtained from the resolution of the non-linear problem, and 
the line follows from our prescription. The inset shows the same quantities in presence of an electrolyte (fCresO = 2.6). 
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FIG. 13. Osmotic pressure (in Pa) as a function of volume fraction. The symbols are the experimental measures of [49] 
for aqueous suspensions of bromopolystyrene particles (with radius a — 51 nm). The continuous curve corresponds to our 
prescription assuming a salt concentration of 10 -6 M in the reservoir. 

In the PB cell model, the osmotic pressure in the solution is related to the densities of micro-ions at the WS cell 
boundary [7,36,47,48]: 

n osm = k B T [p+(R ws ) + p-(Rws) - 2/?o] , (48) 

where we have subtracted the ionic contribution from the reservoir (of salt density p = Jo). This is because the 
electric field vanishes at the WS cell and there is no contribution from the electrostatic pressure at r = Rws- Using 
Eq. (45), Eq. (48) may be recast into 

n osm = ^[^ R (p,/o)-K r 2 es ], (49) 
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where K rcs is the screening constant defined in terms of the ionic strength in the reservoir. Our prescription is supposed 
to give an excellent approximate of the non-linear K pB through ^ LPB , and readily allows an estimation of the osmotic 
pressure. Figure 12 shows the accuracy of the estimate, with or without added electrolyte. 

The comparison of our predictions at saturation to the experimental results reported by Reus et al. [49] is also 
satisfactory, see Fig. 13. It was already pointed out in [49] that PB cell theory reproduced well the experimental data. 
The agreement observed in Fig. 13 however illustrates the relevance of the PB saturation picture -well captured by 
our approach- at high polyion/micro-ion electrostatic coupling (see the discussion in section VII). 

C. Osmotic properties of rod-like polyions 

Expression (49) is also valid in cylindrical geometry, and we show in Fig. 14 the comparison prescription versus 
non-linear PB osmotic pressure. We draw a similar conclusion as for spherical polyions concerning the accuracy of 
our approximation. 




Tl Tl 

FIG. 14. Same as Fig. 12 for cylindrical polyions. Top: salt free suspensions (the inset is a zoom in the small packing 
fraction region). Bottom: situation with added salt (na = 3). 




c c Mol/1 

FIG. 15. Osmotic coefficient of B-DNA solutions as a function of density of DNA phosphate ions c c , for ionic strengths of 
10 mM, 2mM and niM (from bottom to top). The dots are the experimental points obtained from Refs. [50,51], while the 
solid lines correspond to the predictions for <j> using our prescription in cylindrical geometry. The dashed line is the prediction 
of Oosawa-Manning condensation theory. 
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For completeness, we compare in what follows our estimate for the osmotic pressure to the experimental results on 
B-DNA solutions reported in [50,51]. In this work, the authors measured the osmotic coefficient <f> = II osm /II c , defined 
as the ratio between the osmotic pressure II osm to the pressure II C of releasablc counter-ions having bare density c c 
(n c = kBTc c ) against the concentration of B-DNA, a rigid cylindrical polyelectrolyte. A related PB cell analysis may 
be found in [48] while a more thorough investigation has been performed in [52]. 

Within the WS model, B-DNA macro-ions are confined into cylindrical cells, which radius Rws is related to the 
bare concentration of DNA counter-ions as c c — {Idna^Rws)~ 1 -> w ^h &dna = 1.7 A the distance between charges 
along DNA. Note that as in the previous section dealing with charged spheres, there is no adjustable parameter in 
our description since the radius of the DNA and the bare charge (only used to normalize the osmotic pressure to 
II C ) are known from independent measurements. In Fig. 15, the corresponding results for the osmotic coefficient are 
confronted against the experimental data of Refs. [50,51] for various ionic strengths, showing again a good quantitative 
agreement. As in Ref. [7], we report the prediction of classical Oosawa-Manning condensation theory (see e.g. [19,52]), 
for which the osmotic coefficient is constant [<fi = Idna/(^b)] at complete variance with the experiments. In view of 
the results reported in Fig. 14, the disagreement at small c c may be attributed to the (relative) failure of PB theory, 
and not to a weakness of our prescription that should be judged with respect to non-linear PB. 

VII. DISCUSSION - VALIDITY OF THE APPROACH 

Our analysis was carried out at the level of Poisson-Boltzmann theory, which is mean-field like. More refined 
approaches such as the salt free Monte Carlo simulations of Groot [53] for the cell model within the primitive model 
[54] have shown a non-monotonic behaviour of Z c g upon increasing Z for spheres: after the linear regime where 
Z e s ~ Z , Z e g reaches a maximum for a value Z£ aTe and then decreases. When the radius a of the charged spheres 
is much larger than Bjerrum length £b, this maximum is surrounded by a large plateau in excellent agreement with 
the PB saturation value Z sat [15,53]. PB theory appears to become exact for ^b/<i — > [55]. Given that, Z^ are scales 
like (a/^s) 2 and therefore becomes quickly large when the colloid size is increased [53], PB theory is successful in 
the colloidal limit of large a. We recall that this is precisely the limit where our predictions for the effective charge 
at saturation Z sat are reliable (the condition a 3> n^ 1 should a priori be satisfied even if we have shown above that 
our predictions remain fairly accurate down to na of order 1). More generally, the results of Groot [53] show that 
the effective charge Z e g from Monte Carlo simulations within the primitive model for arbitrary a/ is are smaller 
than the quantity Z sat of PB theory. This is a general feature that neglect of ionic correlation (as in PB) leads to 
a underestimated screening of the macro-ion by the micro-ions, and thus to overestimated effective charges (see e.g. 
[52,56] in cylindrical geometry). Our approach thus provides a useful upper limit for a realistic Z c ff. 

A related comment in favor of the validity of PB picture with a saturation plateau for Z c s comes from the work of 
Cornu and Jancovici [57]. For the two dimensional two-component Coulomb gas bounded by a hard wall of surface 
charge a, these authors performed an exact calculation at the inverse reduced temperature e 2 /(kT) = 2 showing that 
the effective surface charge of the wall saturates to a plateau value when a diverges. 

Generally speaking, in an 1:1 electrolyte, PB theory seems to be a reasonable approximation [58,59], all the better 
that the size of the macro- ion is larger than Ib\ the notion of charge renormalization then encaptures the main 
effects of the non-linear PB theory, and is consistent with experimental data on dilute bulk solution [60,61] (see 
also the experimental work cited in section VI). For micro-ions of higher valences (di- or tri-valent), strong ionic 
correlations rule out PB-like approaches, as shown by recent computer simulations of the primitive model [62-64] . As 
a consequence, the results presented here should a priori not be used in the interpretation of experiments involving 
multivalent salts or counter-ions. For a discussion concerning the effects of multivalent counter-ions, we refer to the 
review by Bhuiyan, Vlachy and Outhwaite [65]. 

VIII. CONCLUSION 

The notion of effective charge is widely used in the fields of colloidal suspensions. It allows in practice to describe 
the phase behavior of (highly) charged macro- ions staying at the level of linearized Poisson-Boltzmann equations, 
where the macro-ions are supposed to interact through Yukawa-like pair potentials. However, no general analytical 
description of this renormalization process is available and the effective charge is usually left as a free parameter, 
adjusted to fit experimental (or numerical) data. Physically the charge renormalization process results from the 
strong coupling of the micro-ions in the vicinity of the highly charged colloidal particle. At the level of Poisson- 
Boltzmann theory, the effective charge saturates to a finite value in the limit where the bare charge becomes large. 
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We recall that omission of the non-linearities of PB theory -correctly accounted for by the notion of effective charge- 
may lead to unphysical results (see e.g. [36,66]). 

In the present paper, we have put forward a simple method to estimate the effective charge of highly charged 
colloidal objects either analytically, or through the resolution of a simple equation obtained within linearized Poisson- 
Boltzmann approximation. This approach (mostly suited to describe the colloidal limit ko> 1) amounts to consider 
the highly charged colloids as objects with constant electrostatic potential ~ AkT/e, independently of shape and 
physico-chemical parameters (size, added 1:1 electrolyte...). This result relies on the physical picture that the 
electrostatic energy eVb of the strongly coupled micro-ions (i.e. micro- ions in the vicinity of the highly charged 
macro-ion) does balance their thermal (entropic) energy ksT, resulting in a constant effective surface potential for 
the "dressed" macro- ion. We have successfully tested this approach against a) the geometry of the solid particle, b) 
the confinement (finite concentration situations), c) the presence of added salt, d) exact and approximate solutions of 
the full non-linear PB equations e) direct experimental measurements of the effective charge found in the literature. 
From these different checks, we conclude that our prescription appears to contain the key ingredients involved in 
charge renormalization. 

An important point is that the effective charge is not constant and depends explicitly on the physical conditions 
of the experiment, through ionic strength, density etc .... The effect is quite obvious in the small dilution limit, 
where we found that the (saturated) effective charge is an increasing function of k (for tta > 1, which stems from 
the reduction of the attraction between the counter-ions and the colloid. It pertains for finite concentration and 
the effective charge increases with the ionic strength in the suspension. Addition of salt consequently brings two 
antagonist effects on the effective Coulombic interaction between macro-ions: the range of the interaction decreases 
due to screening, while the amplitude increases due to the effective charge. The competition between these two effects 
might be a key point in the understanding of these systems. It appears therefore interesting to reconsider the phase 
stability of macro- ion suspension in light of these results (see also [68] and more recently [69]). 

Eventually it would be desirable to extend our approach to the case of finite size colloidal particles, such as rods 
with finite length or disks [47] . Accordingly edge effects should show up at the level of our prescription and result in 
an effective charge distribution along the macro-ion, due to the constant potential prescription on the object. Work 
along these lines is in progress. 
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APPENDIX A: ANALYTICAL SOLUTION OF PB EQUATION FOR AN ISOLATED PLANE IN AN 

ELECTROLYTE 



Here, we recall the analytical solution of PB equation for an isolated plane of bare surface charge ae immersed in 
an electrolyte of bulk ionic strength I . In this case, the solution of equation (6) reads [67]: 



— KZ 



where 7 = \J x 1 + 1 — x, k 2 — Stt£b^o and x — k\gc, with the Gouy-Chapman length defined as 

A - = ih- (A2) 

Far from the charged surface, say z > 2k _1 , the solution of PB equation reduces to 

<f> pB {z) * fae-" (A3) 

with (ps = 47. The potential 4*s can be interpreted as the apparent reduced potential extrapolated at contact. In the 
following we shall simply refer to (f>s as the apparent potential. 

As expected, this asymptotic expression for the reduced potential <f> pB in Eq. (A3) precisely matches the solution 
4> LPB of the linearized PB (LPB) equation: 

V Xpb = k Xpb ( A4 ) 

but with the fixed charge boundary condition on the plane replaced by an effective fixed surface potential boundary 
condition </> LPB {z = 0) = 4>s = 47. The effective charge density is then computed using Gauss theorem at the surface, 
yielding 

<7 eff =^. (A5) 

In doing so, we have replaced the initial non linear PB equation with fixed charge boundary condition by the linear 
LPB equation with fixed surface potential. 

Now at fixed k (i.e. constant ionic strength), we progressively increase the bare surface charge a. Accordingly 
kAgc — > and the parameter 7 goes to 1. From Eq. (A5), we obtain that the effective charge and the apparent 
potential <fis have a simple behaviour depending on the comparison of a with cr sat defined as 



cr S at = —r- ■ (A6) 
irl B 



Indeed 



c <C er sat cr e ff = cr 

(j> s ^4<T/<T sat (A7) 
a > c S at a cB = CT sat 

0s = 4 (A8) 

The important point is that in the large bare charge limit, a » (T sat , the effective charge a e s saturates to a value a sat 
independent of the bare one, a. In this limit, the apparent potential also saturate to a constant value, 4>s = 4. 



APPENDIX B: ANALYTICAL SOLUTION OF PB EQUATION FOR A CONFINED PLANE WITHOUT 

ADDED SALT 

An infinite plane (bare surface charge density ae) is placed in the middle of a Wigner-Seitz slab of width 2h. The 
origin of coordinates x = is chosen at the location of the plane such that the volume available to the counter-ions is 
—h < x < h. For symmetry reasons, it is enough to solve the problem for x > 0. The electrostatic potential <j> obeys PB 
equation (4), supplemented with Neumann boundary conditions: V<fi(h) = 0, corresponding to the electroneutrality 
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condition; V</>(0) = — 2itIb(? imposing the charge on the plane. Without loss of generality, we choose the origin of 
potential such that <j>(h) = 0; the analytical solution of PB equation then reads [67] 



i R W = - log 



cos 2 



( {\x\-h) 



1 

PB 



(Bl) 



where K (a) is such that 



hK^„ ( hK m \ 

ZB tan £B = n£ h B2 

V2 V V2 J 

The inverse screening length K pB is related to the density of counter-ions at the WS boundary p~(h) through the 
following expression, reminiscent of the standard definition of the Debye length 

K 2 pB = 4n£ B p-(h). (B3) 

Now we consider the corresponding LPB equation. More precisely, we linearize Eq. (4) around x — h (i.e. the edge 
of the slab). Since we have chosen </> (h) — 0, we impose LPB to vanish at x = h. The resulting equation reads 

V 2 = ^ 2 pB (0 + l) (B4) 

where we have introduced ^ LPB , an "apparent" local Debye screening factor for the linearized PB equation. As 
for the previous PB equation in the no salt case, ^ LPB is not known a priori but results from the electroneutrality 
condition. Indeed, solving Eq. (B4) with the appropriate boundary conditions [V(j)(h) = 0, V0(O) = 47r £b (c/2)] 
yields 

^lpb ( x ) = cosh [^lpb (x-h)]-l, ( B5 ) 

with K LpB (a) such that 

hK LpB sinh [K LpB h] = 2w£ B o-h (B6) 

Note that comparing Eqs. (B2) and (B6), we see that K hpB (a) ^ K pB (a). It is however crucial to remember that 
the LPB solution should not be used with the bare charge a to describe the correct behaviour of (f> pB in the vicinity 
of x = h. 

Next, we implement the procedure proposed by Alexander to find the effective charge in confined situations [5]. 
The effective charge density is accordingly the value of a in the linearized PB equation such that </> pB (x) and </> LPB (x) 
match up to the second derivative at S, the boundary of the WS cell [5]. This condition is equivalent to set 

^ LpB (<7eff) = K pB ( ( 7). (B7) 

Note that in general, whenever the solution of the non-linear PB problem is known, the effective charge cr e ff can be 
directly estimated with Eq. (B7) (this is of course quite academic to obtain in this case an effective charge since 
the full solution for the potential is known; the notion of effective charge is mostly useful in geometries where no 
analytical solution of the PB equation is known). Note also that whenever Eq. (B7) is verified, the third, fourth and 
fith derivative of the linear and non- linear solutions also match at E. 

One deduces eventually the "exact" effective charge from Eq. (B7) for the plane case (by "exact" we mean that 
the effective charge is obtained from the analytical solution of PB equation, in contrast to our prescription) : 

if pB (q) sinh [K pB (a)h] 
CToff= 2^ ' (B8) 

The apparent potential <j)s is also obtained as 

<t>s = LPB (0) = cosh [K pB (a)h\ - 1. (B9) 

From Eq. (B2), we define a critical value for the charge density 
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and we find the asymptotic behaviors 



hK pB = (2a/a c )^ 2 




-3= o> = 5.06er r 



(Bll) 



(B12) 



As in the infinite dilution limit, one obtains that the effective charge a c s coincides with the bare one a for small 
a, and saturates to a finite value when a — > oo. However both a c and the saturation value for the effective charge 
at finite concentration differ from the cr sat of infinite dilution [Eq. (A6)]. We also note that strictly speaking, the 
limits of infinite dilution and of vanishing added salt do not commute: if the limit of vanishing salt is taken first 
(situation investigated in this appendix), before h — > oo, we obtain ^ at = cosh [7r/\/2] — 1, whereas reverting the 



order corresponds to the planar situation of Appendix A with n 
however, the effective charge at saturation vanishes. 



0, and there, we have 



b% at = 4. 



In both cases 
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